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Abstract 

The effects of large CP violating phases arising from the soft SUSY breaking 
parameters on the proton lifetime are investigated in supersymmetric grand uni- 
fied models. It is found that the CP violating phases can reduce as well as enhance 
the proton lifetime depending on the part of the parameter space one is in. Mod- 
ifications of the proton lifetime by as much as a factor of 2 due to the effects of 
the CP violating phases are seen. The largest effects arise for the lightest sparticle 
spectrum in the dressing loop integrals and the effects decrease with the increasing 
scale of the sparticle masses. An analysis of the uncertainties in the determination 
of the proton life time due to uncertainties in the quark masses and in the other 
input data is also given. These results are of import in the precision predictions 
of the proton lifetime in supersymmetric unification both in GUT and in string 
models when the soft SUSY breaking parameters are complex. 



1 Introduction 



It is well known that there are new sources of CP violation in supersymmetric 
theories which arise from the soft SUSY breaking sector of the theory. The normal 
size of such phases is 0(1) and an order of magnitude estimate shows that such 
large phases would lead to a conflict with the current experimental limits on the 
electron and on the neutron electric dipole moment ||]. The conventional ways 
suggested to avoid this conflict is either to assume that the phases are small ||, f§ 
or that the SUSY spectrum is heavy ||. However, recently it was demonstrated || 
that this need not to be the case and indeed there could be consistency with 
experiment even with large CP violating phases and a light spectrum due to an 
internal cancellation mechanism among the various contributions to the EDMs. 
The above possibility has led to a considerable further activity [|7|, [| [J and the 
effects of large CP phases under the cancellation mechanism have been investigated 



in dark matter with the EDM constraints^], in — 2 []TT|1 and in other low energy 
physics phenomena [jT^] . 

In this paper we investigate the effects of large CP violating phases on nucleon 
stability in supersymmetric grand unification with baryon and lepton number vio- 



lating dimension five operators [13, 14, 15]. The main result of this analysis is that 
the dressing loop integrals that enter in the supersymmetric proton decay analysis 
are modified due to the effect of the large CP violating phases. The CP effects 
on the proton life time are most easily exhibited by considering R T defined in 
Eq.(22) which is the ratio of the p lifetime with phases and without phases. R T is 
largely independent of the GUT structure which cancels out in the ratio. Since the 
dressing loop integrals enter in the proton decay lifetime in both GUT and string 
models which contain the baryon and the lepton number violating dimension five 
operators, the phenomena of CP violating effects on the proton lifetime should 
hold for a wide range of models both of GUT and of string variety ]^] . However, 
for concreteness we will consider first the simplest SU(5) supersymmetric grand 
unified model, and then consider a non-minimal extension. As discussed above 
similar analyses should hold for a wider class of models and so what we do below 
should serve as an illustration of the general idea of the effect of large CP phases 
on the proton lifetime. 

The outline of the paper is as follows: In Sec. 2 we give a theoretical analysis of 
the effects of CP violating phases on proton decay in the minimal supersymmetric 
SU(5) model for specificity. In Sec. 3 we discuss the numerical effects of the CP 
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violating phases on R T under the EDM constraints. A non-minimal extension is 
also discussed and an analysis of the uncertainties in the predictions of the proton 
life time due to uncertainties in the quark masses, in (3 P and in the KM matrix 
elements is given. Conclusions are given in Sec. 4. 

2 Theoretical analysis of CP violating phases on 
proton decay in supersymmetric GUTs 

In the minimal supergravity unified model (mSUGRA) flTfl the soft SUSY breaking 
can be parameterized by tuq, mi, Aq, and tan/5, where mo is the universal scalar 
mass, mi is the universal gaugino mass, A is the universal trilinear coupling all 
taken at the GUT scale, and tan/3 =< H 2 > / < Hi > is the ratio of the Higgs 
VEVs where if 2 gives mass to the up quark and Hi gives mass to the down quark 
and the lepton. In addition, the effective theory below the GUT scale contains the 
Higgs mixing parameter which enters in the superpotential in the term [iHxH^- 
In the presence of CP violation one finds that the minimal model contains two 
independent CP violating phases which can be taken to be 6^, which is the phase 
of \i and aA which is the phase of A . 

For more general situations when one allows for non-universalities, the soft 
SUSY breaking sector of the theory brings in more CP violating phases. Thus 
unlike the case of mSUGRA here the U{1) x SU(2) x SU(3) gaugino masses 
(i=l,2,3) can have arbitrary phases, i.e., 

rhi = \mi\e^ (i = 1,2,3) (1) 

While in the universal case a field redefinition can eliminate the common phase of 
the gaugino masses, here one finds that the difference of the gaugino phases does 
persist in the low energy theory and in fact is found to be a useful tool in arranging 
for the cancellation mechansim to work for the satisfaction of the EDMs. In the 
following analysis we carry out an analysis of the proton decay with the most 
general allowed set of CP violating phases. The definition of the mass matrices for 
charginos, neutralinos and for squarks and sleptons have been explicitly exhibited 
in Ref . QTSJ and we refer the reader to this paper for details. The focus of the present 
work is to analyze the effects of CP violating phases on p decay and to estimate 
its size. For the sake of concreteness we begin with a discussion of the simplest 
grand unification model, i.e., the minimal SU(5) model. However, the technique 
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discussed here to include the CP effects on p decay can be used to anlayse the 
CP violating effects for any supersymmetric unified model with baryon and lepton 
number violating dimension five operators. This class includes string models. 

As mentioned above we consider for concreteness and simplicity the minimal 
SU(5) model whose matter interactions are given by JIB], [T4] , |T5fl 



W Y = -lf Uj e uvwxy H?M™Mj y + h, l U- 1 ,,M„M]- (2) 

where M u , M uv are the 5, 10 plet representations of SU(5), and Hi, H 2 are the 5, 5 
of SU(5). After the breakdown of the GUT symmetry and integration over the 
Higgs triplet fields the effective dimension five interactions below the GUT scale 
which governs p decay is given by[[L3|, [L4], [15[ 



Cbl = ^e a UPfiV) l3 (f^^LJ L cAe c Lk (Vu L ) al - u c k d Lal ) + ...) + He. (3) 

Cm = —e^^rhiPVf^ie^u^jUnckdm + •••) + H.c. (4) 

where C^l and C^r are the LLLL and RRRR lepton and baryon number violating 
dimension 5 operators, V is the CKM matrix and fi are related to quark masses, 
and Pi appearing in Eqs.(3) and (4) are the generational phases given by = (e l7i ), 
£i7i = (i=l,2,3). 

The baryon and the lepton number violating dimension five operators must be 
dressed by the chargino, the gluino and the neutralino exchanges to generate ef- 
fective baryon and lepton number violating dimension six operators at low energy 
(Some examples of dressing loop diagrams are given in Fig.l). It is in this process 
of dressing of the dimension five operators that the CP violating phases of the soft 
SUSY breaking sector enter in the proton decay amplitude. The CP phases enter 
the dressings in two ways, via the mass matrices of the charginos, the neutralinos 
and the sfermions, and via the interaction vertices. Taking account of this addi- 
tional complexity, the analysis for computing the proton decay amplitudes follows 
the usual procedure. Thus to dress the dimension five operators the squark and 
slepton fields must be eliminated in terms of their sources. As an example, the 
up squarks in the presence of CP violating phases can be eliminated using the 
relations 

u lL = 2 J [A L ul L ui + Af R R m ] 

u iR = 2j [A*R ui + Af L L m ] (5) 
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where 

T - 5Ll R - 5Ll 

Here L/ is the sum of fermion-sfermion-gluino, fermion-sfermion-chargino and 
fermion-sfermion-neutralino interactions and 

= [\D uin \ 2 A uil + \D uil2 \ 2 Au i2 ) 
A* = [\D ul21 \ 2 A uil + \D ui22 \ 2 Au i2 \ (7) 

and 

A^ = —D u mD U ii2[A U ii — Aui 2 ] 

A* L = -D um D* uil2 [A uil - AtZfi,]. (8) 

Here Un and are the squark mass eigenstates for the squark flavors Un and Ui 2 
and A Ui \ and A Ui2 are the corresponding propagators, and D ui is the diagonalizing 
matrix for the Ui squarks, i.e., 

DlM 2 p u% = diagiM 2 ^ M 2 J (9) 

We note the special arrangement of the complex quantities and their complex 
conjugates in Eqs.7 and 8. Specifically we note that while in the absence of CP 
phases A^ = A^ L this is not the case in the presence of CP phases and in general 
one has A^ ^ A^ L as is seen from Eqs.(8). L ui and R ui defined by Eq.(6) receive 
contributions from the chargino, the neutralino and the gluino exchanges. 

Following the standard procedure fi~3|, [14] , |15| one obtains the effective dimension 
six operators for the baryon and the lepton number violating interaction arising 
from dressing of the dimension five operators. From this effective interaction one 
obtains the proton lifetime decay widths for various modes using the effective 
Lagrangian methods. We limit ourselves here to the dominant decay mode p — > 
ViK + . Including the CP violating effects the decay width for this process is given 
by 



P 2 m N 



-> 



M 2 H ^fl 



m 



N 



;i + mjv( 3 ^ 3F) )(i + yf + (e-^y~ 9 + w 2 + §y?s*) 

+1— d(i + yf - (e-^y § - y 2 )5 i2 + ^l^s)! 2 (10) 

6 TUB Ao 
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where 

A ViK = (sin2/3M^)- 1 a 2 2 P 2 m c m^ i t i^i^ 22 [^(q d t ; W) + F(c; e<; W)] (11) 

In the above Al(As) are the long (short) suppression factors, D,F, f n are the effec- 
tive Lagrangian parameters, and (3 P is defined by (3 p Ul = t a b c t a p < 0\d2 L u^ L u2 L \p > 
where is the proton wavefunction. Theoretical determinations of (3 P lie in the 
range 0.003 — 0.03 GeV 3 . Perhaps the more reliable estimate is from lattice gauge 
calculations which gives |T9[1 j3 p = (5.6 ± 0.5) x 10~ 3 GeV 3 . 



Aside from the explicit CP phases via the exponential factor e - ^ 3 in Eq.(10), 
CP effects enter dominantly in .F's which are the dressing loop integrals. For the 
chargino exchange in the presence of CP violating phases one has 

J A=l,2 

G A [A^U A1 -A L /ep* A2 ) (12) 

Here G A (A=l,2) are the propagators for the chargino mass eigenstates and the 
matrices U and S enter in the biunitary transformations to diagonlize the chargino 
mass matrix Mq such that 

U*M C S^ = diag(m x +,m x +) (13) 

In Eq.(10) the quantities y\ k are the corrections due to the chargino exchanges 
involving third generation squarks, 3^ is the contribution from the gluino exchange, 
y% is the contribution from the neutralino exchange, and are the contributions 
from the dressing of the RRRR dimension 5 operators. 
The gluino exchange contribution y§ is given by 

V - 4 -Pi a 3 rn u V n H(u; d; g) - H(d; d; g) , , 

~ 9 ~ 3 P 2 a 2 m c V 21 V 2 \V 22 F(c; s; W) + T{c; faW) 

The function TC is defined by H(u;d;g) = f(rh u ; rhd, rh g ) where f is defined by 
Eq.(19) below. The contributions y^ from the dresssing of the RRRR dimension 
five operators are given by 

R = Pi mtmiVnVvVL Q(f;i;W) 

1 P 2 m c m h V 21 V22V^ x T{c 1 b-W)+T{c-f-W) 



and by 



R = Pi m t m s V^V 12 V^ Q(f;t;W) 

2 P2m c m b V 2l V22V^J r {c~b-W)+J r {c-f-W) 
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where Q's are defined as follows: 
Q(f;t;W0 = -327A( ™ T - ) / £ A^ 2 [e t Af ^ 2 - Af S* A1 ]G A (17) 

\/ ZMu/ COS 1 J a -, o 



A=l,2 



The dressing loop integrals can be expressed in terms of the basic integral 

f(/jLi, fij, Hk) = -16tt 2 z J AiAjGk (18) 
where 2222 

f '{f^ii [^j 1 ^k) 2" 2~[ 2 2^ 2 2" 2^ 2"] (^) 

— ^fe A*i — f^j ~ ^k ^k 

In the limit of no CP violation the analysis limits correctly to the previous results 
which do not include CP violating effects. 

3 Numerical effects of CP violating phases on 
proton decay 

We discuss now the numerical size of the effects of CP phases on the proton 
lifetime. As is obvious from our discussion above the proton lifetime is highly 
model dependent. Specifically there are two main factors that govern the lifetime 
of the proton. One of these depends on the nature of the GUT sector, i.e., if the 
GUT group is SU(5), SO(10), E 6 ,..etc and on the nature of the GUT interaction, 
e.g., on the GUT Higgs structure, while the second factor that controls proton 
decay is the sparticle spectrum and the sparticle interactions that enter in the 
dressing loop integrals. 

In the simplest SU(5) GUT model with two Higgs multiplets Hi and H 2 , GUT 
physics enters mainly via the Higgs triplet mass M# 3 . The proton decay lifetime is 
significantly affected if we were to change the GUT structure. Thus, for example, 
if one had many Higgs triplets |20|, Hi and Hi (i=l,..,n) where only the triplets Hi 



and Hi couple with matter, i.e., 

W triplet = g iJ+ J Hi + HiMijHj (20) 

then the effective interaction on eliminating the Higgs triplets is 

Wt ff = -J(M n V (21) 

Here one finds that the effective Higgs triplet mass is = ((M _1 )n) _1 . The 

model with M# 3 and would have very similar CP violating effects for the 
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same low energy sparticle spectrum. The reason is that the CP effects are largely 
governed by the nature of the low energy physics, e.g., the sparticle mass spectrum 
and the couplings of the sparticles with matter. Thus we expect similar size CP 
violating effects in models with different GUT structures but with similar size 
sparticle spectrum 

To discuss the CP violating effects on the proton lifetime it is useful to consider 
the ratio R T (p —* v + K + ) defined by 

Rrip ^v + K + ) = t( P ^v + K + )/r (p ^u + K + ) (22) 

Here r(p — > v + K + ) is the proton lifetime with CP violating phases and r (p — > 
u + K + ) is the lifetime without CP phases. This ratio is largly model independent. 
Thus most of the model dependent features such as the nature of the GUT or 
the string model would be contained mostly in the front factors such as the Higgs 
triplet mass, the quark masses, the A s and A L suppression factors all of which 
cancel out in the ratio. Similarly the quantity (3 P which is poorly known cancels 
out in the ratio as do the KM matrix elements. 

We analyze R T under the constraints that CP violating phases obey the exper- 
imental limits on the electron and the neutron EDMs. For the electron and for 
the neutron the current experimental limits arejl], |2| 

|4| < 4.3 x l(T 27 ecm, \d n \ < 6.3 x 1(T 26 ecm (23) 

We are interested in the effects of large phases on the proton lifetime and for 
these to satisfy the EMD constraints we use the cancellation mechanism. In Fig. 2 
we present five cases where for different inputs the electron EDM is plotted as a 
function of 9^. An analysis of the neutron EDM for the same input is given in 
Fig. 3. One finds the cancellation mechanism produces several regions where the 
EDM constraints are satisfied. In Fig. 4 we give a plot of R T for the same set of 
inputs as in Figs. 2 and 3. The analysis shows that R T is a sensitive function of the 
CP phase 9^ and variations of a factor of around 2 can occur. We also note that 
both a suppression as well as an enhancement of the proton lifetime can occur as a 
consequence of the CP violating effects. Interestingly the largest CP effects on R T 
occur here at the points of maximum cancellation in the EDMs as may be seen by a 
comparison of Figs. 2, 3 and 4. The variations in R T due to the phases arise because 
of constructive and destructive interference between the exchange contributions of 
chargino 1 (xt) an d chargino 2 (xt) ( see Fig.l). We give an illustration of this 
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phenomenon in Table 1. The analysis of Table 1 exhibits the cancellation in the 
imaginary part of the amplitude for the decay process p — > v + K + from chargino 
1 and chargino 2, and this cancellation leads to an enhancement in the p lifetime 
ratio for this case. 

It is possible to promote each of the cancellation points in Figs. 2 and 3 into a 
trajectory in the m — m.i plane by scaling upwards by a common scale transformation] 

m — > Am , mi — > Ami (24) 

2 2 

The size of the sparticle spectrum depends on the scale A. In general, the larger 
the value of A the heavier is the sparticle spectrum and correspondingly smaller 
is the CP effect on the dressing loop as demonstrated in Fig. 5. For the minimal 
SU(5) case one needs a relatively heavy spectrum with some of the sparticle masses 
~ 1 TeV to stabilize the proton which has the current experimental limit for the 
p — > vK decay mode of r(p — > uK) > 5.5 x 10 32 yr[21]. Because of the heaviness of 



the sparticle spectrum, the CP effects for the minimal SU(5) model are typically 
small, i.e., of order only a few percent. Larger CP effects can occur in non-minimal 
models where one has several Higgs triplets. Thus we consider an example where 
one has two pairs of heavy Higgs triplets with the Higgs triplet mass matrix given 
by 



/ aA 
\aA M 2 



(25) 



Such a structure can arise, for example, in an SO(10) model|22| with two 10's 
of Higgs and a 45 of Higgs with a superpotential of the type Wh = M 2 102 H + 
A10ih45 h 10 2 h- After the 45 of Higgs develops a VEV < 45 H >=(a, a, a, 0, 0) x 
%a 2 one finds that only one pair of Higgs doublets remain massless while the Higgs 
triplets (H t i,H t2 ) and (H t i,H t2 ) have the mass matrix given by Eq.(25). In this 
case one has M e // = a 2 A 2 /M 2 and one can arrange for proton stability even with 
a light spectrum by an adjustment of the parameters aA and M 2 . 

Finally we discuss the current uncertainties in the proton lifetime predictions. 
Uncertainties arise from the errors in the quark masses, in (3 P and in the KM 
matrix elements. The largest source of uncertainties arises from the strange quark 
mass (m s ). There are several determinations of m s : m s = 193 ± 59 MeVP3], 
m s = 200 ± 70 MeVEl, m s = 170 ± 50 MeVp51,m s = 155 ± 15 MeVpq, all 



evaluated at the scale lGeV. We take for our average m s = 180 ±50 MeV. For the 
charm quark mass (m c ) we use m c = 1.4 ±0.2 GeV^J while for the bottom quark 
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mass (m&) we use m& = 4.74 ± 0.14 GeV||27|| . The contributions from the first 
generation quarks are small and are not the sources of any significant uncertainty 
in the p lifetime. The errors in the KM matrix elements are of a subleading order 
for the vK mode but are still significant enough to be included. We use the results 
of Ref.[2~7 for the allowed ranges of the KM matrix elements. For (3 P we use the 



result of the lattice gauge analysis of Ref. . In Fig. 6 we exhibit the error corridor 
for the proton lifetime for the case (1) of Fig.2 with M 2 /aA = 0.01, M 2 = Mq- 
One finds that given the current errors in the input data the predictions for the 
proton lifetime has an uncertainty of about a factor of 2 (1^05) on either side of 
the mean. A similar analysis holds for the cases(2-5) of Fig.2. We note that the 
uncertainties in the predictions of the proton lifetime is of the same order as the 
size of the CP violating effects. It is for this reason that we choose to exhibit 
the results of our analysis in Figs. 4 and 5 in terms of the ratio R T since the 
effects of the uncertainties cancel in the ratio. The analysis also shows that an 
improvement in the determination of the quark masses and of (5 P is essential for 
a more precise prediction of the proton lifetime in supersymmetric unification of 
the type discussed here. The reduction of the error in the prediction of p lifetime 
will also help to define the CP effects on proton decay when such a decay is 
experimentally observed. 

In summary the CP violating effects on the proton lifetime are relatively large if 
the sparticle spectrum entering the dressing loop integrals is relatively light and the 
CP violating effects get progressively smaller as the scale of the sparticle spectrum 
entering the dressing loops gets progressively larger. The current experimental 
limits on the sparticle masses allow for a relatively light sparticle spectrum, i.e., 
significantly smaller than 1 TeV. This means that there exists the possibility of 
significant CP violating effects on the proton lifetime. However, the minimal 
SU(5) model does not support the scenario with a light spectrum and thus the 
CP violating effects for the case of the minimal model are small. However, for the 
non-minimal case proton stability can occur even for a relatively light spectrum 
due to suppression from a more complicated Higgs triplet sector. In these types 
of models CP violating effects can be significant. 

4 Conclusion 

In this paper we have investigated the effects of CP violating phases arising from 
the soft SUSY breaking sector of the theory on the proton decay amplitudes. It is 
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found that the CP effects can increase or decrease the proton decay rates and that 
the size of their effect depends sensitively on the region of the parameter space one 
is in. Effects as large as a factor of 2 are seen to arise from CP violating phases 
in the part of the parameter space investigated and even larger effects in the full 
parameter space may occur. It is found that the CP violating effects in the min- 
imal SU(5) model are typically small since a relatively heavy sparticle spectrum 
is needed to stabilize the proton in this case and a heavy spectrum suppresses the 
CP effects in the dressing loop integral. However, significantly larger CP effects 
on the p lifetime are possible in non-minimal models with more than one pair of 
Higgs triplets since in these models the proton can be stabilized with a relatively 
light sparticle spectrum. We also investigated the uncertainties in the p lifetime 
predictions due to uncertainties in the quark masses, in (3 P and in the KM matrix 
elements. We find that these uncertainties modify the proton lifetime by a factor 
of 2 around the mean value. The observations arrived at in this analysis would 
be applicable to a wide class of models, including GUT models and string models 
with dimension five baryon and lepton number violating operators 



Note Added: After the paper was submitted for publication an improved limit on 
p — > v^K + mode of r(p — ► u fl K + ) > 1.9 x 10 33 yr has been reported||28||. The new 
limit does not affect the conclusions arrived at in this paper. 
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Table 1:CP effects on Chargino dressings. 



case 


Chargino 1 


Chargino 2 


(Re<7IAJ, Im^/IAJ) 


(-.22, -.89) 


(-.026, .2) 


case 


(ReA v J\A v J, \mA u J\A u J) 


\AuJ/\A„J 


sum 1&2 


(-.25, -.69) 


.74 



Table caption: The table gives an analysis of the dressing loop integrals for 
dresssings with Charginos 1&2 (see Fig.l) and their sum for the case when 
m = 71 GeV, mi = 148 GeV, tan/? = 2, 9^ = 1.4, & = 0.3, f 2 = 1-8, & = 0, 
where all phases are in radians. The analysis shows cancellations in the dressings 
between Chargino 1 and Chargino 2 for the case with phases. 



Figure Captions 

Fig.l: Examples of the dressing of LLLL baryon- number- violating dimension five 
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operators by chargino, gluino and neutralino exchanges that contribute to the pro- 
ton decay. Cancellation among diagrams such as between xl an d xt exchanges 
can lead to an enhancement of the proton lifetime. The dressings of the RRRR 
dimension five operators is also exhibited. 

Fig. 2: Plot of Logio|d e | vs 9^ exhibiting cancellations where the five curves cor- 
respond to the five sets of input for the parameters tan/5, mo, mi/2, £i, £2, £3, 
a Ao , and A given by (1)2,71,148, -1.15,-1.4,1.27, -.4,4 (dotted), (2)2,71,148, 
-.87,-1.0,1.78, -.4,4 (solid), (3)4,550,88, .5,-1.55,1.5, .6, .8 (dashed), (4)4,750,88, 
1.5,1.6,1.7, .6, .8 (long dashed), and (5)2,71,148, .55,1. ,1.35, -.4,4 (dot-dashed). All 
masses are in GeV and all phases are in radians. 

Fig. 3: Plot of Logio|<i„| vs 9 ^ exhibiting cancellations where the five curves corre- 
spond to the five sets of input for the parameters tan/5, m , m^, £1, £2, £3, ®a , 
and A as given in Fig. 2. 

Fig. 4: The ratio R T of the proton lifetime with phases and without phases as a 
function of 9^ for the five cases given in Fig.2. 

Fig. 5: The ratio R T as a function of the scaling factor A defined in the text. The 
four curves correspond to the four sets of input for the parameters tan f3, £1, £2; £,3, 
9^, o:a () and A given by (1)2,— 1.15,— 1.4,1.27,— 1.7,— .4,4 with m = 71 and mi/ 2 = 
148 for the point of intersection with R T axis (dotted)., (2)2, -.87,-1. 0,1. 78, -2. 15, -.4,4 
with m = 71 and mi/ 2 = 148 for the first point (solid). (3)4, .5, — 1.55, 1.5, 1.56, .6, .8 
with m = 550 and mi/ 2 = 88 for the first point (dashed). (4)4,1.5,1.6,1.7,— 1.56, .6, .8 
with m = 750 and mi/ 2 = 88 for the first point (long dashed). All masses are in 
GeV and all phases are in radians. All trajectories satisfy edms constraints. 
Fig. 6: Exhibition of the uncertainties in the proton lifetime predictions due to 
uncertainties in the input data for case(l) of Fig.2 where we assumed M 2 /aA = 
0.01, M 2 = Mq. 
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